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Abstract 

Spin correlations of top quarks produced in hadron collisions have not been observed experimen- 
tally with large significance. In this Letter, we propose a new variable that may enable demonstra- 
tion of the existence of spin correlations with 3 — 4 a significance using just a few hundred dilepton 
events both at the Tevatron and the LHC. Such number of dilepton events has been observed at 
the Tevatron. At the LHC, it will become available once integrated luminosity of a few hundred 
inverse picobarns is collected. 



The existence of spin correlations of top and 
anti-top quarks in tt pair production in hadron 
collisions is a solid prediction of the Standard 
Model. The possibility to observe these cor- 
relations is unique to top quarks since their 
large masses, short lifetimes and the relative 
weakness of chromomagnetic fields in the QCD 
vacuum, make it difficult for non-perturbative 
effects to depolarize t and t before they de- 
cay. Therefore, if top quarks are produced in a 
particular polarization state, spin correlations 
can be observed by studying kinematic distri- 
butions of the top quark decay products which 
are sensitive to t and i polarizations. For ex- 
ample, in the dilepton channel pp(pp) —> ti —> 
bb 11 vv, the V — A structure of the charged cur- 
rent forces momenta of anti-leptons (leptons) 
to be aligned (anti-aligned) with the direction 
of the top (anti-top) spin vectors. 

The traditional way to study top quark spin 
correlations 0, |, i, i, 1, 1, 0, 1, i, M ^ _ f airry 
complex. It involves choosing the t and t spin 
quantization axes and identifying suitable ref- 
erence frames and angular distributions that 
are sensitive to these correlations. Because of 
the unobserved neutrinos in the dilepton events, 
full kinematics can not be reconstructed and 
determination of quantization axes and refer- 
ence frames becomes difficult. This feature and 
a relatively low yield of dilepton events at the 
Tevatron is partially responsible for the fact 
that spin correlation measurements performed 
by the CDF and DO collaborations are not con- 
clusive. For example, the parameter n related 



to the top quark spin asymmetry in the dilep- 
ton channel at the Tevatron is predicted with a 
very small uncertainty in the Standard Model, 
k = 0.78 P, 0, S I, |, 0, S 0, \M- How- 
ever, it is measured to be n — O.llg^l and 



K = 0.3±° ;| by DO [111 and CDF [lj collabo- 
rations, respectively, with 5.4(3.0) fb -1 of inte- 
grated luminosity. Although these results are 
consistent with the Standard Model, they do 
not demonstrate the existence of top quark spin 
correlations with sufficient significance. A sim- 
ilar situation occurs when spin correlations are 
measured in the lepton plus jets channel [13 1. 

It is then natural to ask if a better way ex- 
ists to establish the presence of spin correla- 
tions convincingly. This question was recently 
discussed by G. Mahlon and S. Parke in Ref. 
[14] ]. They suggested that spin correlations at 
the LHC can be observed by measuring the 
relative azimuthal angle Acf> of the two lep- 
tons from top decays in the laboratory frame, 
provided that only events with the low invari- 
ant mass of tt pairs, M t i < 400 GeV, are ac- 
cepted. While in this case it is possible to dis- 
tinguish spin-correlated and spin-uncorrelated 
events, Ref. Q recognizes that placing a cut 
on M t t is unphysical since, in dilepton events, 
the tt invariant mass can not be fully recon 



structed on an event-by-event basis. Ref. [14 1 
suggests that one can put a cut on the statisti- 
cally reconstructed invariant mass M t j but this 
cut does not seem to work as well as the cut on 
Mu proper. It was later shown in Refs. [HEil 
that simpler cuts on kinematics of final state 
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particles - for example an upper cut on the 
transverse momentum of the charged leptons - 
lead to the A<j> laboratory frame distributions 
that are sufficiently different to enable distin- 
guishing between spin-correlations and no-spin- 
correlations hypotheses. 

While Ref. [14j opened up a new direction in 
the studies of top quark spin correlations, simi- 
lar toprevious papers on the subject [H, [2, H, 0, 
HiESfiSEll it focused on a single kinematic 
distribution. However, since many kinematic 
features of a particular event may be sensitive 
to top quark spin correlations, we should try 
to use all the information present in a particu- 
lar event to establish their existence. To that 
end, we ask if, given a set of ti events observed 
at the Tevatron or the LHC, it is possible to 
distinguish the hypothesis that spins of the ti 
pair, entangled in the production process, re- 
main entangled at the time of their decay, from 
the hypothesis that strong QCD dynamics de- 
polarizes produced top quarks and kinematic 
features of their decay products are not corre- 
lated. 

It is well-understood by now that the op- 
timal way to answer this question requires the 
construction of a likelihood function, where the 
optimality should be understood in the sense 
of Neyman-Pearson lemma 17]. In our case, 
this lemma explains how to minimize the prob- 
ability to accept the spin-correlation hypothesis 
when it is false, for fixed probability to reject 
spin-correlation hypothesis when it is true. We 
discuss how to construct the likelihood func- 
tion. However, before going into this, we point 
out that we do not consider issues of exper- 
imental resolution in this Letter. In reality, 
experimental resolution may be important, so 
our results should be considered as the estimate 
of the highest significance with which two hy- 
potheses can be separated. 

We study dilepton events that contain two 
neutrinos in the final state, pp(pp) — > ti — > 
bb vi^v^. Because neutrino momenta can 
not be measured, the number of measurable 
kinematic variables £c b s is smaller than the 
number of kinematic variables x = (a; bs, ^unobs) 
needed to describe a particular event. The prob- 
ability distribution for events with £c b s set to a 
particular value is computed by integrating dif- 
ferential cross-sections for a hypothesis H over 



unobserved kinematic variables. We write 

P H (x ohs ) = Mh 1 J ckr unobs aH d ^ X> , (1) 

where Afn is the normalization factor. Given 
the two hypotheses, that the top quarks spins 
are correlated (H = c) or uncorrelated (H = 
u), we introduce a variable, related to a likeli- 
hood ratio for a single event, that emphasizes 
the difference between the two hypotheses 

7Z(x b-) — Pc{ x obs) ^ 

-fc(^obs) ~\~ -^ti(^obs) 

We then calculate the probability distribution 
of the likelihood variable 7?.(a; bs), given a par- 
ticular hypothesis about the underlying physics 

p H (K) = a£ J dx^^5(K(x) - K), (3) 

and perform statistical tests to see how many 
events are required to achieve the separation of 
the two hypotheses H — c,u. 

It is important to realize that, at the ex- 
pense of claiming that our likelihood ratio 1Z 
is the optimal observable |18| to separate spin- 
correlation and no-spin-correlation hypotheses, 
we can use different cross-sections to construct 
the likelihood variable lZ(x) in Eas. (|ll2p and 
to calculate the probability distribution pu (7Z) 
in Eq.([3]). We note that we use the Born dif- 
ferential cross-section da^(x) to define 7?.(x). 

This is a good choice because da^ (x) captures 
main kinematic features of the actual physi- 
cal process and it is inexpensive computation- 
ally. However, since this choice does not corre- 
spond to the actual probability distribution of 
the dilepton events, strictly speaking, TZ is not 
the optimal variable. Nevertheless, as long as 
1Z helps to separate the two hypotheses, opti- 
mality is not essential. We emphasize, however, 
that we use the best available approximation to 
the true cross-section day (a:) to construct the 
realistic probability distribution of the variable 
1Z. To this end, in this Letter we employ the 
next-to-leading order (NLO) QCD prediction 
for the top pair production do~ t i that includes 
top quark spin correlations and radiative cor- 
rections to top quark decays [l9j . 

Since we use the leading order cross-section 
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Figure 1: Normalized distributions of the likeli- 
hood variable R for the spin-correlation and no-spin- 
correlation hypotheses at the 7 TeV LHC. The NLO 
QCD results for the distributions are shown. The 
bands correspond to the choices of the renormaliza- 
tion/factorization scales 0.5 rat < fJ. < 2mt- 



Figure 2: Normalized distributions of the likeli- 
hood variable R for the spin-correlation and no-spin- 
corrclation hypotheses at the Tevatron. The NLO 
QCD results for the distributions are shown. The 
bands correspond to the choices of the renormaliza- 
tion/factorization scales 0.5 nit < M < 2mt. 



to compute 1Z, the following issue appears. In 
general, the NLO QCD approximation includes 
processes with additional massless particles in 
the final state. Therefore, we need a prescrip- 
tion of how to map the kinematic features of 
such final states onto leading order kinematics. 
Indeed, at leading order the process pp(pp) — > 
ti — > bb 11 viv\ has two massless 6-quarks. As- 
sociating these 6-quarks with two &-jets recon- 
structed according to a well-defined jet algo- 
rithm solves the problem of additional radia- 
tion in the event. However, perturbatively, b- 
jets at leading order are massless, while this 
is not necessarily true in higher orders. This 
feature makes it difficult to connect the lead- 
ing order kinematics that enters the calculation 
of 1Z with kinematics of the actual event. To 
address this problem, we adopt the Ellis-Soper 
jet algorithm [2(j, where reconstructed jets are 
always massless. 

The discussion in the previous paragraph 
tells us how to map kinematics of a higher- 
order process to the kinematics of a tree-level 
process. As input for the calculation of the 
likelihood TZ(x h s ), we use four-momenta of the 
two &-jets, the four-momenta of the two charged 
leptons and the missing transverse momentum, 
which we identify with the component of the 
momentum of the two neutrinos, orthogonal 
to the collision axis. We also note that, since 
charges of &-jets can not be unambiguously de- 
fined, we require a procedure to assign one jet 
to be a 6-quark jet and the other jet to be a 



6-quark jet. We do this by computing the in- 
variant mass of the positively charged lepton 
and the two 6-jets and identifying the jet that 
minimizes this invariant mass, with the &-quark 
jet. The other 6-jet is then identified with the 
&-quark jet and, in leading order kinematics, we 
treat this jet as if it comes from the decay of 
the anti-top quark. 

Having discussed a procedure to identify the 
input, we turn to the calculation of 1Z; this re- 
quires an integration of the tree-level differen- 
tial cross-section for pp(pp) — > ti — > bb 11 viUi 
over unobserved components of the neutrino 
momentum. In general this is difficult, but we 
assume that the process goes through the on- 
shell intermediate states, so that the invariant 
masses of blv and blv are equal to m t and that 
the invariant masses of Iv and lv are equal to 
my/. Hence, we compute 

PH{Pobs,P±Aniss) = N H 1 Jl&PvW&Pv] 

obs j Pv j Pv ) | 
ij (4) 

x 8{Ml - m 2 w )5{Mf, - m 2 w ) 
x5(Ml b ~mj)S(M^ b -m 2 t ), 

where [dp] = d 3 p/((27r) 3 2_B) is the invariant 
integration measure, p obs = {Pb,Ph,Pi,Pi} is 
the set of observable momenta, fi(x) are par- 
ton distribution functions and k = (pi + 



pj + .. +Pk) 2 are the respective invariant masses 
squared. As we see, there are six (5-functions 
in Eq.Q, so that all integration variables are 
fixed; all we need to do is to solve the on- 
shell constraints. This is a standard procedure 
which is described e.g. in Ref. [2~lj : we do not 
repeat such a discussion here. In general, solv- 
ing the on-shell constraints leads to several so- 
lutions (the maximal number is four), in which 
case all these solutions should be taken into ac- 
count. 

The result for the probability distribution 
reads 

-Pff (Pobs,RL,miss) = N H X y^y^ J a 

ij a (5) 

where the second sum is over all the solutions 
that are obtained by reconstructing the final 
state and J a is the Jacobian which appears 
when the integration over the neutrino momen- 
tum is carried out. Also, /jj = fi{x^\) is a 
parton distribution whose argument is recon- 
structed from the kinematics of a final state. 
Finally, we emphasize that in the calculation 
of probability distributions Ph and the vari- 
able TZ, we always use leading order matrix el- 
ements, as explicitly shown in Eq.©. 

The result for Ph in Eq. ([5]) allows us to cal- 
culate the likelihood TZ and carry out the indi- 
cated program. In practice, however, we make 
use of the fact that both at the Tevatron and 
the LHC there is a single partonic channel that 
dominates the production process. Therefore, 
in Eq.([S]), we use i = (u, d), j = (u, d) to com- 
pute TZ for the Tevatron and i = j = g to com- 
pute TZ for the LHC. We also neglect the depen- 
dence of the normalization factor Mh on the 
hypothesis H, following the observation that 
total cross-sections for tt pair production are 
insensitive to the (non)existence of top quark 
spin correlations. Finally, in the computation 
of the likelihood variable TZ{x \, s ), we always 
set the renormalization and factorization scale 
to m t and use leading order parton distribution 
functions. 

To calculate the probability distribution of 
the variable TZ for a given hypothesis H, we 
perform a numerical integration in Eq.Q. We 
generate events assuming that they must pass 
basic selection cuts for the tt events. For both 
the LHC and the Tevatron, we require p l ± > 
20 GeV, pf ss > 40 GeV and \rji\ < 2.5. We 



also require that there are at least two 6-jets in 
the event; jets are defined using Ellis-Soper jet 
algorithm [H with AR = \J An 2 + A<t> 2 = 0.4 
and the jet p^-cut is set to 25 GeV for both the 
Tevatron and the LHC. We use MRST2001 and 
MRST2004 parton distribution functions [H 
in LO and NLO computations, respectively. All 
calculations that we report in this paper make 
use of the numerical program for computing 
NLO QCD effects in tt pair production, devel- 
oped in Ref. [11]. 

We now present the results of the calcula- 
tion. First, we compute the distribution of the 
likelihood variable TZ for both spin-correlations 
and for no-spin-correlations hypotheses. These 
distributions are shown in Figs lll2l for the LHC 
and the Tevatron, respectively. The two distri- 
butions are similar although the LHC distri- 
bution is more narrow. Also, as follows from 
Figs lll2l the scale dependence of ^-distributions 
is small and we neglect it in what follows. 

For both the Tevatron and the LHC, there 
is a difference between the two distributions, 
which is especially visible in the region of small 
TZ. To find the number of events that is re- 
quired to distinguish between the two TZ distri- 
butions, we perform a statistical test [23(. To 
this end, we generate N events according to 
the probability distribution ph(TZ) defined in 
Eq.Q and calculate the quantity 

L = 2 ln[£ c /£ u ], (6) 

JV 

where Ck = Yi Pk(Ri)- The statistical inter- 
pretation of L can be found in Ref. [23[ • We re- 
peat this procedure multiple times, for H = c,u 
and obtain two distributions of the variable L. 
The distribution of the variable L is peaked 
at positive (negative) values if events are gen- 
erated with the hypothesis H = c (H = u) 
since, on average, C c > C u (C u > C c ). Exam- 
ples of such distributions are shown in Fias l3l4l 
To compute the significance S with which the 
hypotheses H = c and H = u can be sepa- 
rated, we find the point beyond which the right- 
side tail of the left histogram and the left-side 
tail of the right histogram have equal areas. 
These areas correspond to the one-sided Gaus- 
sian probability outside of the S/2a range. If 
the two L-distributions are Gaussian with unit 
widths, the significance S is the separation be- 
tween peaks of the two distributions. 

The significance with which two hypothe- 
ses can be separated depends on the number 
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Figure 3: Distributions of the likelihood ratio for cor- 
rect/false hypotheses, assuming five hundred dilepton 
events at the 7 TeV LHC. Arbitrary units. We show 
LO and NLO QCD results for the L-distribution com- 
puted with fjj = rat- 



Figure 4: Distributions of the likelihood ratio for cor- 
rect/false hypotheses, assuming three hundred dilepton 
events at the Tevatron. Arbitrary units. We show LO 
and NLO QCD predictions for the L-distribution com- 
puted with fi = rat- 



al events N with which the two hypotheses are 
probed. To understand what is a reasonable 
value of N, we note that the pp —> ti produc- 
tion cross-section at the -J]i = 7 TeV LHC is 
approximately 160 pb 24j,[25[. Since W^-bosons 
decay to electrons and muons twenty percent 
of the time, and assuming thirty percent effi- 
ciency, we find that 1 fb _1 of the integrated 
luminosity corresponds, roughly, to 2500 dilep- 
ton events. It is expected that 1 fb" 1 of lumi- 
nosity will be collected at the LHC by the end 
of 2011 and this sets reasonable upper bound 
on the number of leptons N. 

In fact, we do not need that many. We take 
N = 500, which corresponds to 200 pb" 1 , as- 
suming 30% efficiency. We then consider 10 6 
pseudo-experiments and obtain the two distri- 
butions shown in FigJ31 We convert the over- 
lap of the two distributions into statistical sig- 
nificance and find that, with 500 events, the 
two distributions shown in Fig[TJ can be sepa- 
rated at the 4<t level. It is interesting to note 
that the difference between NLO and LO L- 
distributions at the LHC is very small, cf. FigEl 

We now turn to the discussion of the ti 
production at the Tevatron. The production 
cross-section of the it pairs at the Tevatron is, 
approximately, 7 pb (for the latest measure- 



ments, see Refs. [26J, l27j). Taking the accu- 



mulated luminosity to be 6 fb _ and assum- 
ing 30% efficiency, we find that five hundred 
dilepton (/i, e) events at the Tevatron should 
have been observed. We take N — 300 and, 



by considering 10 6 pseudoexperiments, we ob- 
tain L-distributions shown in FigSl I n this 
case, there are significant differences between 
L-distributions computed at leading and next- 
to-leading order. Analyzing the L-distribution 
obtained with the NLO QCD approximation, 
we find that, with 300 Tevatron dilepton events, 
the spin-correlation hypothesis can be estab- 
lished with the significance that is close to 3.5<r. 

Summary: We have shown that a likelihood- 
based analysis should make it possible to demon- 
strate the existence of top quark spin correla- 
tions in dilepton events at the Tevatron and the 
LHC. We constructed the relevant likelihood 
function and computed its probability distri- 
bution through next-to-leading order in per- 
turbative QCD. Neglecting all the experimen- 
tal uncertainties and the background contribu- 
tions that are relatively small for the dilepton 
channel, we find that with 500 dilepton events 
at the LHC and with 300 dilepton events at 
the Tevatron the existence of spin correlations 
can be established with better than 3<r signifi- 
cance. This number of events will require just 
about 200 pb -1 accumulated luminosity at the 
LHC and is already available at the Tevatron. 
We believe that our results are sufficiently en- 
couraging to warrant a more complete study 
including proper treatment of experimental un- 
certainties and backgrounds. 
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